Tunneling-percolation origin of nonuniversality: theory and experiments 
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A vast class of disordered conductirig-insulating compounds close to the percolation threshold 
is characterized by nonuniversal values of transport critical exponent t, in disagreement with the 
standard theory of percolation which predicts t ~ 2.0 for all three dimensional systems. Various 
models have been proposed in order to explain the origin of such universality breakdown. Among 
them, the tunneling-percolation model calls into play tunneling processes between conducting par- 
ticles which, under some general circumstances, could lead to transport exponents dependent of the 
mean tunneling distance a. The validity of such theory could be tested by changing the parameter a 
by means of an applied mechanical strain. We have applied this idea to universal and nonuniversal 
Ru02-glass composites. We show that when t > 2 the measured piezoresistive response V, i. e. the 
relative change of resistivity under applied strain, diverges logarithmically at the percolation thresh- 
old, while for t ~ 2, F does not show an appreciable dependence upon the Ru02 volume fraction. 
These results are consistent with a mean tunneling dependence of the nonuniversal transport expo- 
nent as predicted by the tunneling-percolation model. The experimental results are compared with 
analytical and numerical calculations on a random-resistor network model of tunneling-percolation. 

PACS numbers: 72.20.Fr, 64.60.Fr, 72.60.+g 



I. INTRODUCTION 

Despite that transport properties of disordered 
insulator-conductor composites have been studied for 
more than thirty years, some phenomena still remain in- 
completely understood. One such phenomenon is cer- 
tainly the origin of nonuniversality of the DC transport 
near the conductor-insulator critical transition. Accord- 
ing to the standard theory of transport in isotropic perco- 
lating materials, the bulk conductivity a of a composite 
with volume concentration x of the conducting phase be- 
haves as a power-law of the formic 

a ~ o- (x - x c y, (1) 

where ao is a proportionality constant, x c is the critical 
concentration below which the composite has zero con- 
ductivity (or more precisely the conductivity of the insu- 
lating phase) and t is the DC transport critical exponent. 
The above expression holds true in the critical region 
x — x c <C 1 in which critical fluctuations extend over dis- 
tances much larger than the characteristic size of the con- 
stituents. As a consequence, contrary to Co and x c which 
depend on microscopic details such as the microstructure 
and the mean inter-grain junction conductance, the ex- 
ponent t is expected to be material independent i£ The 
universality of t is indeed confirmed by various numeri- 
cal calculations of random-resistor network models which 
have established that t = to oi 2.0 for three dimensional 
lattices to a rather high accuracyi^^ 

Confirmations to the standard percolation theory of 
transport universality are found only in a limited num- 
ber of experiments on real disordered composites. This is 
illustrated in Fig^ where we have collected 99 different 
values of the critical exponent t and the critical thresh- 
old x c measured in various composites including carbon- 
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FIG. 1: Collection of critical exponent values t and cor- 
responding critical threshold concentration x c for various 
disordered conductor-insulator composites. Carbon-black- 
polymer systems are from Refs.|al2ql. oxide-based thick film 
resistors are from Refs. |27M38j . metal- inorganic and -organic 

insulator granular metals are from Refs. mmm. The 

dashed line denotes the universal value to ~ 2.0. 

black-polymer systems^ - — oxide-based thick film re- 
sistors (TFRs)j2£~— and other metal- inorganic and - 
organic insulator composites^^ - — It is clear that, de- 
spite that many of the i-values reported in Fig[I]are close 
to to ~ 2.0, almost 50% of the measured critical expo- 
nents deviate from universality by displaying t ^ to. 
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Examining all the data reported in Fig^ one observes 
that the lack of universality is not limited to a particular 
class of materials, although the granular metals (empty 
diamonds in Fig^l have somewhat less spread values of 
t compared to the carbon-black and TFRs composites. 
Another important observation is that for the vast ma- 
jority of the cases, the nonuniversal critical exponent is 
larger than and only few data display t < to. Finally, 
there is no clear correlation between t and the critical 
concentration value x c . 

With the accumulation of experimental reports of 
nonuniversality, various theories have been proposed in 
order to find an origin to this phenomenonAiS^fliSiiS In 
Ref@ it was argued that, in carbon-black-polymer com- 
posites, long-range interactions could drive the system 
towards the mean-field regime for which t = < m f = 3.0. 
This interpretation cannot however explain the observa- 
tion of critical exponents much larger than t m { , such as 
those of carbon-based composites or TFRs which display 
values of t as high as t ~ 5 — 10. The authors of Refl5Cl 
introduced the random-void (RV) model of continuum 
percolation where current flows through a conducting 
medium embedding insulating spheres placed at random. 
By using an earlier result^ they were able to show that 
for this model DC transport is described by an univer- 
sality class different from that of standard percolation 
model on a lattice. The resulting critical exponent was 
found to be t ~ 2.4 for three dimensional systems. More 
recently, Balberg generalized the RV model in an attempt 
to explain higher t- values^ The same author proposed 
also a model of transport nonuniversality based on an 
inverted RV model in which current flows through tun- 
neling processes between conducting spheres immersed in 
an insulating medium^ Within this picture, if the distri- 
bution function of the tunneling distances decays much 
slower than the tunneling decay, then the critical expo- 
nent becomes dependent of the mean tunneling distance 
a and in principle has no upper bound. 

In addition to the above models, there were also expli- 
cations pointing out that when Eq. Q is used to fit exper- 
imental data not restricted to the critical region, "appar- 
ent critical exponents", usually larger than to, could be 
misinterpreted as real critical exponentsi^LS Although 
this possibility cannot be excluded for some of the data 
reported in Fig^ it is however quite unrealistic to iden- 
tify the whole set of reported nonuniversal exponents as 
merely apparent. 

The mean-field interpretation^ the RV model and its 
generalization's^*^ and the tunneling inverted RV model 
(also known as the tunneling-percolation model), 51 have 
been devised to describe nonuniversality for various 
classes of materials. For example, the RV model applies 
in principle to composites where the linear size of the 
conducting particles is much smaller than that of the in- 
sulating grains so that the conducting phase can be ap- 
proximated by a continuum. The tunneling-percolation 
model has been instead conceived to apply to those com- 
posites for which inter-grain tunneling is the main micro- 



scopic mechanism of transport, such as in carbon-black- 
polymer composites, 54 or in oxide-based TRFsi 5 ^ 5 ^ 5 ^ In 
principle, the two models could even coexist together if 
the continuum phase of a RV system is made of non- 
sintered conducting particles interacting through tunnel- 
ing processes. 

In this situation, the different proposed theories could 
account for the variety of nonuniversal exponents shown 
in FigQ] However, it is also true that, in order to identify 
a given mechanism of nonuniversality for specific com- 
posites, little has been done beyond a mere fit to Eq.QJ. 
For example, In Refl52l the DC critical exponents have 
been examined together with the relative resistance noise 
exponent, and in Refl39l a study on AC and magnetore- 
sistive exponents for various composites has been pre- 
sented. As a matter of fact, no conclusive answers have 
been reached and the different models of nonuniversal- 
ity listed above have not proven to really apply to real 
composites. 

In this paper we present our contribution to the un- 
derstanding of the origin of transport nonuniversality by 
attacking the problem from a different point of view. 
In contrast with the mean-field hypothesis^ the RV 
models and its extension^ the tunneling-percolation 
model of Balberg predicts that the critical exponent t 
acquires an explicit dependence upon a microscopic vari- 
able (the mean-tunneling distance a) which could be al- 
tered by a suitable external perturbation. So, if transport 
nonuniversality is driven by tunneling, it would be possi- 
ble to change the value of the transport critical exponent 
t by applying a pressure or a strain to the composite. 
Conversely, when a material belongs to some universality 
class (standard percolation theory, mean-field universal- 
ity class, or the RV model) its exponent is expected to be 
independent of microscopic details and an applied strain 
would not change t. 

We have applied this idea to Ru02-based TFRs whose 
transport properties are known to be governed by inter- 
grain tunneling processes ££££££ In the following of this 
paper we show that the behavior of the piezoresistive 
response, i. e. the change of resistivity upon applied 
mechanical strain^ as a function of concentration x of 
Ru02, can be interpreted as due to a tunneling dis- 
tance dependence of the DC critical exponent i, as orig- 
inally proposed in Refl5ll The paper is organized as fol- 
lows. In the next section we briefly review the tunneling- 
percolation theory and the RV model and its extension. 
In Sec lIIII we describe the theory of piezoresistivity for 
percolating composites and in Sec II VI we present our ex- 
perimental results. The last section is devoted to a dis- 
cussion and to the conclusions. 



II. MODELS OF NONUNIVERSALITY 

The RV modelsj£2i&£ and the tunneling-percolation 
theory^i have one point in common. They all rely on the 
work of Kogut and Straley who first developed a theoret- 
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ical model of nonuniversality based on random-resistor 
networks^ In this model, to each neighboring couple of 
sites on a regular lattice it is assigned with probability 
p a bond with conductance g ^ and bond with g = 
with probability 1 — p. The resulting bond conductance 
distribution function is then: 



p(9)=ph(g) + (l-p)S(g), 



(2) 



where 5(g) is the Dirac delta function and h(g) is the dis- 
tribution function of the finite bond conductances. Close 
to the bond percolation threshold p c , the conductivity £ 
of the network behaves as: 



^o{p-PcY 



(3) 



where S is a prefactor. In this and in the subsequent 
section we distinguish the conductivity S of a random- 
resistor network from that of real composites [Eq.flJ]. 
When h(g) = 5(g — go) where go is some nonzero value, 
Eq. (J2J reduces to the standard bi-modal model for which, 
close to the percolation threshold, the network conduc- 
tivity £ follows Eq.@ with t = t ~ 2.0 for all three 
dimensional lattices. Instead, if h(g) has a power law 
divergence for small g of the form: 



lim h(g) oc g' 



(4) 



where a < 1, then universality is lost for sufficiently large 
values of the exponent a~ For lattices of dimension D, 
the resulting conductivity critical exponent ig^SiSLSil 



t 



t if (D-2)u+j^ < t Q 

(D-2) V+J ^ if (D-2) V+J ^>to 



(5) 



where to is the universal value and v is the correlation- 
length exponent (v — 4/3 for D — 2 and v ~ 0.88 for 
L> = 3). For D — 3 and by using t ~ 2.0 and v ~ 0.88 
the value of a beyond which universility is lost is a c — 
1 - l/(i — y) — 0.107. Only recently Eq.© has been 
demonstrated to be valid to all orders of a e = 6 — D 
expansion in a renormalization group analysis^ 

In the original work of Ref 49], a was considered as 
no more than a parameter of the theory without a justi- 
fication on microscopic basis. This came later with the 
RV models and the tunneling-percolation theoriesi^SiiS 
However, before discussing the microscopic aspect, we 
find it interesting to point out that Eq. (|5J) predicts that t 
cannot be lower than to and it is in principle not bounded 
above. This is in qualitative agreement with the experi- 
mental values oft reported in Fig^ Furthermore, Eq.JSJ) 
can give us some information on the distribution of t val- 
ues expected by the theory. In fact, given some normal- 
ized distribution function f(a) for the parameter a and 
for D = 3, the distribution N(t) of the t values is: 



N(t) 



daf(a)6[t - t(a)] = S(t - t ) / daf(a) 



+ / daf(a)6 



t - v 



1-a 



(6) 
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FIG. 2: Distribution of the t-values reported in Fig^ The 
solid line is Eq.(|7J renormalized in order to reproduce the 
total number of data. Inset: log-log plot of the distribution 



with a fit to the power law a(t 
b = 2.5 ± 0.4. 



with a = 192 ± 62 and 



If we assume that f(a) can be approximated by a con- 
stant / for a c < a < 1, then the above expression reduces 
to: 



N(t) = 11- 



I 



5(t-t Q )+f 



1 



0(t-to), (7) 



which predicts a rapid decay N(t) oc (t — v)~ 2 for t > t - 
Equation JJJ) is plotted in Fig. (J2J together with the distri- 
bution of t values reported in Fig.JQ). The distribution 
N(t) has been renormalized to the number of t values 
and the constant / has been fixed to reproduce the num- 
ber of data with t > 3.0. There is an overall qualitative 
agreement between the distribution of the experimental t 
values and N(t). In particular, the asymmetry of the dis- 
tribution and its tail for t > to — 2.0 are well reproduced. 
Fitting to a power law leads to a decay proportional to 
(t — J /)- 2 - 5 ±°- 4 ( S ee inset of Fig|3J| which is in fair agree- 
ment with the predicted behavior (i — v)~ 2 . Due to the 
limited number of t values available to us, the agreement 
with Eq.© could be fortuitous. However, the point here 
is that the gradual decrease of the number of times high 
values of t are reported is not necessarily due to "bad 
fits" to Eq.QfSLS but it can be, at least qualitatively, 
explained by the form of Eq. © • 

Let us now discuss the microscopic origin of the ex- 
ponent a. According to the original RV model^S the 
crucial parameter is the conducting channel width 5 left 
over from neighboring insulating spheres. For D = 3, the 
cross-section of the conducting channel has roughly the 
shape of a triangle and the resulting channel conductance 
g scales as g oc <5 3 / 2 . 50 Hence, if p(S) is the distribution 
function of the channel width <5, the distribution h(g) of 
the conductances reduces to: 



h(g) = I dS P (S)S(g - 50 <S 3/2 ) 



(8) 
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where go is a proportionality constant. For random dis- 
tribution of equally sized spheres, p(S) is a constant for 
S — ► and Eq.JHJ) gives fr(<?) oc g -1 / 3 for g — > Oi^ /i(#) is 
therefore of the same form of Eq.QJ with a = 1/3 > a c 
which, according to Eq.JSJ, predicts a critical exponent 
t = v + | ~ 2.38. Despite that this value is only slightly 
larger than the universal exponent to ~ 2.0, the RV 
model has the merit of being a fully microscopic justi- 
fication of Eq. 10} . In order to allow for higher values of 
t, in ReflH^ it has been proposed to relax the condition 
that p(S) is a constant for S — ► by introducing the more 
general condition p oc 8~ u where u> < 1. By using Eq.lJEJ 
one readily finds that h(g) is again given by Eq.0J but 
j£2m Hence, the critical exponent is t = to 



with a = | 
for u> < lo c 



>«<■ 



+ |/(1 - u) for 



, | ~ 0.4 or t = v + | y 
w > oj c , i. e., t is no bounded above. The use of p oc 5~ u 
has been justified to be a simple ansatz to describe real 
composites in which correlations between the conducting 
and insulating phases may yield to deviations from the 
ideal RV system. 

It is important to stress that the RV model of Refl50l 
does not give a breakdown of universality, but rather to 
a universality class different from that of the standard 
DC transport percolation on a lattice. In particular, the 
exponent t — does not depend on the conductivity of 
the continuum and an applied isotropic strain or pressure, 
albeit affecting the channel width 6, does not change the 
relation g oc <5 3 / 2 . This should remain true also for the 
generalized RV theory of ReflH^. 

Contrary to the RV models, the tunneling-percolation 
theory of nonuniversality allows for a transport critical 
exponent dependent of the microscopic conductivites. 51 
In this model, current flows through tunneling processes 
between neighboring conducting spheres dispersed in an 
insulating medium and, for sufficiently low concentra- 
tions of the conducting spheres, the ensemble of tunnel- 
ing bonds form a percolating network. The coexistence 
of tunneling and percolation has been recently settled by 
experiments probing the local electrical connectivity of 
various disordered systems, and a recent review on this 
issue can be found in Ref|6j In what follows, we con- 
sider the situation in which grain charging effects can be 
neglected with respect to the tunneling processes, as en- 
countered in systems with sufficiently large conducting 
grains and/or high temperatures. Let us consider then 
the inter-particle tunneling conductance: 



g = g e { > /t; , 



(9) 



where go is a constant, £ is the tunneling factor which 
is of the order of few nm, r is the distance between the 
centers of two neighboring spheres of diameter $ (r > $ 
for impenetrable particles). If we denote with P(r) the 
distribution function of adjacent inter-sphere distances r, 
the bond conductances are then distributed according to: 



Kg) 



drP{r) 5 



g- 9oe 



-2(r-*)/S 



Balberg notices that if P{r) had a slower decay than 
Eq.@ for t — ► oo, then h(g) would develop a divergence 
for g — > 0. For example, by assuming that 



P(r) 



-(r-*)/(tt-#) 

a - $ 



(11) 



is a good approximation for the (normalized) distribution 
of inter-particle distances, then Eci. (|10J) would reduce to 



Kg) 



go 



a I g 

go 



(12) 



with a = 1 — where a is the mean distance between 

a— <P 

<r i n 

— — ^ 1 — ™ ~ 0.9, transport 



2 ^ M. It is impor- 



(10) 



neighboring particles. For < 1 — a r 
becomes nonuniversal with t = v 

tant to stress that here nonuniversality is not driven by 
geometrical factors as in the RV model, but rather by 
physical parameters such as £ and a. These can be dif- 
ferent depending on the composite and can be modified 
by a suitable external perturbation. In fact, in the case 
of a applied pressure or strain, the mean tunneling dis- 
tance a would change leading to a modification of the 
critical exponent t. As already pointed out in the intro- 
duction, the detection of such an effect would be a direct 
signature of a tunneling-percolation- like mechanism of 
transport nonuniversality. 

The tunneling-percolation theory in its original for- 
mulation relies on Ea. (|llf) which should be regarded as 
a phenomenological model of the distribution function 
of inter-particle distances. Recently, however, a micro- 
scopic derivation of Eci. (|llf) has been formulated for a 
bond-percolation regular network in which the bonds 
have probability p of being occupied by a string of n 
non-overlapping spheres^ In this case, in fact, Ea. l|ll|) 
is the exact distribution function for spheres placed in 
a one-dimensional channel^ and the resulting bond 
conductance distribution is proportional to g~ a ' 1 where 
a n = 1 — n _ $ and a n is the mean inter-particle dis- 
tance for a ffond occupied by n spheres. Despite of its 
over-simplification, this construction shows however that 
the tunneling-percolation mechanism of nonuniversality 
can be justified by a fully defined model, without need of 
phenomenological forms of P(r). 



III. THEORY OF PIEZORESISTIVITY 

The sensitivity of the tunneling-percolation mechanism 
of nonuniversality to variations of the mean tunneling dis- 
tance can be exploited by imposing a volume compres- 
sion or expansion to the system. As we discuss below, 
under these circumstances the relative change of resis- 
tivity, i. e. the piezoresistive response, changes dramati- 
cally depending on whether the DC exponent is universal 
(t = to) or instead it is driven away from universality by 
the tunneling-percolation mechanism. 
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Let us consider the rather general situation in which 
a parallelepiped with dimensions L x , L y and L z is sub- 
jected to a deformation along its main axes x, y, and 
z. The initial volume V = L x L y L z changes to V(l + </>) 



where <f> — e x + e y + e z is the volume dilatation and 
Si = 5Li/Lt are the principal strains along with 
i — x,y, z. In the absence of strain, we assume that 
the conductivity £ of the parallelepiped is isotropic, so 
that the conductance Gi measured along the i axis is 
Gi = YiLjL^I Li. For small e, =/= (i = x,y,z), the 
conductance variation SGi is therefore: 



SGi 
Gi 



£ 



where 



5£i 
£ 



r±(e 3 -+e fc ) 



(13) 



(14) 



are the relative variation of the conductivity along the 
i — x,y, z directions. The coefficients T|| and T± are the 
longitudinal and transverse piezoresistive factors defined 
as: 



dln(Ei) _ d\n(pi) 



r, = - 



de l 
dei 



dei 

dHPj) 
dei 



(15) 

(i * 3), (16) 



where pi — £^ is the resistivity along the i-axis and In is 
the natural logarithm. The distinction between longitu- 
dinal, T||, and transverse, T±, piezoresistive responses is 
important whenever the values of the strains Si depend 
on the direction, as it is encountered when the sample 
is subjected to uniaxial distortions as those induced in 
cantilever beam experiments (see next section). 

In what follows, we are mainly interested in the 
isotropic (or hydrostatic) piezoresistive factor T defined 
as the resistivity change induced by the isotropic strain 
field £i = e for all i = x,y, z: 



r = - 



dln(S) d\n{p) 



de 



de 



(17) 



which can be obtained by applying an hydrostatic pres- 
sure to the parallelepiped. However, T can also be ob- 
tained by setting Si — e in Eg. 1)14(1 yielding: 



r = r 



2r 



(18) 



which is a useful relation when the experimental set up 
does not permit to apply an isotropic strain field. 

Let us now study how the tunneling-percolation the- 
ory of non-universality affect the piezoresistive factor T. 
To this end, we assume that a cubic bond-percolation 
network is embedded in a homogeneous elastic medium 
and that the elastic coefficients of the network and the 
medium are equal. Under an isotropic strain field £j = e 
(i = x, y, z) the mean tunneling distance a changes to 
a(l + e) independently of the bond orientation. Hence, 



by assuming for simplicity that <& — > <I>(1 + e), the tun- 
neling parameter a = 1 — -^i entering Eq.JSJ) becomes 
a —> a + (1 — a)e for e <C 1. According to the discus- 
sion of Sec ITTland to Eas. (|3l5ll7|) . close to the percolation 
threshold p c , the piezoresistive factor behaves therefore 
as: 

a < a c 

• *]n(p-p c ) a>a c ' (19) 



where 



To 

dt 
de 



rfln(£ ) dln(£ ) 



de 



da 



I- a 



= 2(a -*)/£, 



(20) 
(21) 



where we have used da/de = (1 — a). The tunneling dis- 
tance dependence of the DC transport exponent is there- 
fore reflected in a logarithmic divergence of T as p — > p c . 
Instead, when a < a c , the DC exponent remains equal 
to to also when e ^ and the resulting piezoresistive fac- 
tor is simply equal to T — To, independently of the bond 
probability p. 

It is worth to note that, as shown in Refl64l contrary 
to r, the breakdown of universality has no effect on the 
piezoresistive anisotropy defined as \ = (rii — r_i_)/r|| 
which behaves as x oc (p — p c ) x where the critical expo- 
nent A is independent of a also when a > a c ~ 

Equation (|19fl is an exact result as long as we are con- 
cerned with the p — p c dependence close to the perco- 
lation threshold. However, in addition to the prefactor 
of the logarithm, T depends also on the tunneling pa- 
rameter a through the term To- This dependence is far 
to be trivial. In fact, consider a tensile strain (e > 0) 
which enhances the bond tunneling resistances leading 
to an overall enhancement of the sample resistivity. In 
this case, T — d\n(p)/de must be strictly positive. This 
means that To > when a < a c while, from Eci. (|19f) . 
when a > a c To does not need to be positive to en- 
sure r > because of the presence of the logarithmic 
divergence. In the next subsections we provide evidences 
that indeed To changes sign in passing from a < a c to 
a > a c by showing that both the effective medium ap- 
proximation and numerical calculations on cubic lattices 
give negative values of To for nonuniversal DC transport. 
Together with the logarithmic divergence of T for p —> p c , 
a negative value of To would be an additional signature 
of a tunneling-percolation mechanism of nonunivcrsality. 



A. effective medium approximation 

In the effective medium approximation (EMA) the con- 
ductivity £ = g/i of a bond percolation cubic lattice 
(with bond length £) is obtained by the solution of the 
following integral equation (see for example Ref 2): 



dgp(g) 5 — g = 0, 
3 + 2ff 



(22) 
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FIG. 3: Prefactor Eo of Eq.lJ^J as a function of the tunnel- 
ing parameter a. Note that Eo is a decreasing (increasing) 
function of a for a < (a > 0). Inset: the p independent 
contribution To to the piezoresistive factor P. Note that Eo 
and To are not plotted for small values of a because in this re- 
gion the logarithmic corrections to E calculated within EMA 
change the simple power law behavior of Eq.©. 



where p(g) is the distribution function of the bond con- 
ductances g given in Eq.J2J. Close to the percolation 
threshold p c (p c — 1/3 in EMA) and by using the 
tunneling-percolation distribution of Eq. II 21) with go — 1 , 
the above expression reduces to: 



(1 - a)g / dg 



9 + 2.g 



Pc)- 



(23) 



It is clear that for g — ► the integral in Ea. (|23ll remains 
finite as long as a < 0, while for a > it diverges as 
g~ a . Hence, for p — > p c the conductivity £ follows the 
power law behavior of Eq.© with t = 1 for a < and 
t = 1/(1 — a) for a > 0. The corresponding prefactor So 
can be evaluated explicitly: 



S o (a<0) = 



S (a > 0) 



3 a 



21 -a 



*7(a)7(2 — a) 



(24) 



(25) 



where 7 is the Euler gamma function. Note that the 
above expressions gives S (a — > 0) = 0, which is a result 
due to enforcing £ to behave as Eq.©. Actually, at 
a = and for g < 1 Eq.lO reduces to —gln(2g) = 
p c ) which leads to logarithmic corrections in the p — p c 
dependence of g. For a^0, the logarithmic corrections 
are not important only in a region around p = p c which 
shrinks to zero as a — > 0. 

S (a) is plotted in Fig|3|as a function of a. For a < 
(a > 0), So( a ) is a decreasing (increasing) function of 
a. Hence, according to the second equality of Ea.p0)l. 
To is expected to be positive for a < and negative for 
a > 0. This is confirmed in the inset of FigEI where 



Tq is plotted as a function of a. Note that as a — > 1, 
To goes to —00. In fact, from Eq. l|25ll . in this regime 
Sq (a) ~ ±3 1 / ( - 1 - a '> = |3* which implies that: 



r n 



-ln(3) 



efe 



(26) 



B. Monte Carlo calculations on cubic lattices 

To evaluate the a dependence of the prefactor So of 
Eq. we have used the transfer-matrix method applied 
to a cubic lattice on N — 1 sites in the z direction, N sites 
along y, and C along the x direction^ Periodic bound- 
ary conditions are used in the y direction while a unitary 
voltage is applied to the top plane, and the bottom plane 
is grounded to zero. For sufficiently large C (£ 3> Af) 
this method permits to calculate the conductivity Sjv 
per unit length of a cubic lattice of linear size TV. The 
transfer matrix algorithm is particularly efficient at the 
bond percolation threshold p c ~ 0.2488126, and it is usu- 
ally used in connection with finite size scaling analysis of 
S^v to extract highly accurate values of the DC exponent 

a 

In performing the calculations we have considered the 
following geometries: = 6 (£ = 5x 10 7 ), N = 8 (C — 
2 x 10 7 ), N = 10 (£ = 1 x 10 7 ), N = 12 (£ = 8 x 10 6 ), 
N = 14 (£ = 4 x 10 6 ), and N = 16 (£ = 2 x 10 6 ). 
From Eq.©, the resulting conductivity S^ for finite N 
at p = p c can be written as: 



S JV = S (a)[p c -p c (iV)] 



t(«) 



(27) 



where we have explicitly written the a dependence of the 
prefactor and of the exponent. In the above expression, 
Pc{N) is the percolation threshold of a finite systems of 
linear size N. Only at A^ — > 00, p c {N) coincides with 
the percolation threshold p c of an infinite system, while 
for finite values of A^ the two quantities are related via a 
finite size scaling relation of the typeiSi 



Pc - p c (N) ~ AN-V U (1 + BN~ U ), 



(28) 



where v ~ 0.88 is the correlation length exponent, A and 
B are constant and u is the first scaling correction ex- 
ponent. All the constants appearing in Ea. (|28fl depends 
solely on the connectivity of the network, and are there- 
fore independent of the tunneling factor a. By inserting 
Ed. 1(2*8)1 in Eq.J37|), the conductivity reduces to: 



where 



: N ~ ^ (a)N-^ a ^ v (l + BN- W )^ a \ (29) 



S (a) = So(a)A t («). (30) 



Our strategy to calculate S (a) is the following. We 
first fit or numerical data of Sjv with Eq. I|29|) by setting 
uo fixed. In this way we obtain the exponent t{pt) and 
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FIG. 4: DC transport exponent t calculated from the transfer- 
matrix method as a function of the parameter a. The solid 
line is the exact result Eq.©. 
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FIG. 5: Prefactor So (a) of the finite size scaling relation 
Eg. 1291 as a function of a. Inset: behavior of Eo(a) as a func- 
tion of log( — a)/[l + log(— q)] (where log is the logarithm to 
base 10) for a = -1, -10, -100, -1000, and a = -co. This 
latter case corresponds to a Dirac delta distribution function: 
Hd) = for whi ch we obtain E (-oo) = 0.960 ±0.007. 



the prefactor So (a). We repeat this procedure for var- 
ious values of a ranging from a ~ 1 down to a = — oo 
which corresponds to the Dirac delta distribution func- 
tion h(— oo) = 8(g — 1). In this limit, So (—00) is known 
with a rather good accuracy, permitting us to calculate 
from Eq. 1|30|) the value of the constant A. In this way we 
finally obtain So(a) = So(a)/A*( Q ' by using the values 
of the exponent calculated before. 

In FigQJwe plot the values of the exponent t as a func- 
tion of a obtained by setting to = 1 in Eq.JSHJ). We have 
checked that this choice for u produced the best overall 
agreement with the exact result Eq.Q shown in FigQ] 
by the solid line. In accord with the a-independence 
of B. we noticed little deviations from B ~ 0.7 in the 
whole range of a. The agreement between the calcu- 



FIG. 6: Prefactor Eo(a) of the conductivity obtained from 
Eq.JSOl with A = 1.55 and t from Fig^J. Inset: semi- 
logarithmic plot of Eq (a) as a function of the DC exponent 
t. Solid line is a fit with E (q) = afe* with a = 0.018 ± 0.004 
and b = 1.9 ±0.1. 



lated exponent and Eq.JSJ) is very good far away from 
a = a c ~ 0.107. In the vicinity of a c the competition 
between two different fixed points leads to a less good 
agreement, as already noticed in previous works^i 

In Fig|3]we report the calculated values of the prefactor 
So (a) of the finite size scaling relation Ea. H29l) . Note that 
the overall dependence of So (a) upon a resembles that 
of Fig|3 although the presence of A^ ' certainly affects 
the a > a c region. In the inset of Fig[S] we have plot- 
ted the behavior of So (a) as a — > —00. In this regime, 
the exponent is universal (t = to ~ 2) and the whole 
cv-dependence is contained in the conductivity prefactor 
So(a). From Ref 5 we know that So(— 00) ~ 0.4, while 
we have obtained So(— 00) — 0.96. Hence from Ea. l30|) 
we have A ~ (0.96/0.4) 1 / 2 ~ 1.55. 

Our final results for S (a) = S (a)/1.55^ Q ^ are plot- 
ted in FigEI where the t(a) values are those plotted in 
FigEl As for the EMA case, So(cv) decreases for a suf- 
ficiently smaller than a c ~ 0.107 while it increases for 
a > a c . Hence, also our Monte Carlo calculations con- 
firm that the p-independent part To of the piezoresistive 
response is positive or negative, depending whether a is 
less than or larger than cv c , respectively. In the inset 
of Fig|S]we report a semi-logarithmic plot of So (a) as a 
function of the DC exponent t. For high values of t, the 
data are reasonably well fitted by a relation of the form 
S (a) = a&* with a = 0.018 ± 0.004 and b = 1.9 ± 0.1 
(solid line in the inset of FigHjJ). Hence: 



r = -]n(fc)^~-0.6ffc 
de de 



(31) 



confirming the asymptotic formula Ea. H26|) obtained 
within EMA. 
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IV. EXPERIMENT 

In this section we describe our experiments aimed 
to investigate the piezoresistive response of disor- 
dered conductor-insulator composites made of conduct- 
ing Ru02 particles embedded in a insulating glass. As 
shown in FigQ] this kind of TFRs displays both universal 
(t ~ to) and nonuniversal (t > to) behaviors of transport, 
although the factors responsible for such changes have 
not yet been identified. In addition, transport in such 
kind of materials is known to be governed by electron 
tunneling through the glassy film separating two neigh- 
boring conducting particles £££§£1 Hence, Ru02-based 
TFRs are ideal systems to test whether a tunneling- 
percolation mechanism of transport nonuniversality sets 
in. 

Our samples were prepared starting with a glass frit 
with the following composition: PbO (75% wt), B2O3 
(10% wt), SiC>2 (15% wt). In order to avoid crystallisa- 
tion, 2% wt of A1 2 3 was added to the glass powder. Af- 
ter milling, the glass powder presented an average grain 
size of about 3 /zm, as measured from laser diffraction 
analysis. Thermogravimetric measurements showed neg- 
ligible loss in weight, indicating glass stability and no 
PbO evaporation during firing up to 800 °C, and differen- 
tial scanning calorimetry measurements indicated a soft- 
ening temperature of about 460 °C. For the conductive 





FIG. 7: SEM images of the surface of the A series with RUO2 
volume fraction x = 0.08 and nominal RUO2 grain size of 
400 nm (see text) for different firing temperatures Tf. (a): 
Tf = 525° C; (b):7> = 600° C 
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FIG. 8: SEM images of the surface of the B series with Ru02 
volume fraction x = 0.08 and nominal RUO2 grain size of 
40 nm (see text) for different firing temperatures Tf. (a): 
T s = 550° C; (b):Tf = 600° C 



phase we used two different Ru02 powders with nom- 
inal grain sizes of 400 nm (series A)and 40 nm (series 
B). Transmission electron microscope analysis confirmed 
that the finer powder was made of nearly spherical par- 
ticles with a diameter of about 40 nm, while the coarser 
powder had more dispersed grain sizes (100 nm mean) 
with less regular shape. 

TFRs were then prepared by mixing several weight 
fractions of the two series of Ru02 powder with the glass 
particles. An organic vehicle made of terpineol and ethyl 
cellulose was added in a quantity of about 30% of the 
total weight of the Ru02-glass mixture. The so-obtained 
pastes were screen printed on 96% alumina substrates on 
pre-fired gold terminations. For the conductance mea- 
surements, eight resistors 1.5 mm wide and different 
lengths ranging from 0.3 mm up to 5 mm were printed 
on the same substrate. The resistors were then treated 
with a thermal cycle consisting of a drying phase (10 min 
at 150 °C) followed by a plateau, reached at a rate of 20 
°C/min, of 15 min at various firing temperatures Tf (see 
belowand TableQJ. After firing, the thickness of the films 
was about 10 /im. 

In Fig[7| and Fig|S| we show scanning electron micro- 
scope (SEM) images of the surfaces of A and B series, 
respectively, with Ru02 volume concentration x = 0.08. 
The firing temperature T f is T f = 525°C and Tf = 600°C 
for the A series [Figs. Ufa) and (b), respectively], while 
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Tf = 550°C and Tf = 600°C for the B series [Figs. Ufa) 
and (b), respectively]. In the images, dark areas indicate 
highly conducting regions rich of RuC>2 clusters, while 
the white zones are instead made of insulating glass. The 
grey regions surrounding the RuC>2 clusters indicate that 
some conduction is present, although much lower than 
the dark areas, which we ascribe to finer Ru02 particles 
dispersed in the glassiS At the length scale shown in 
the figures the conducting and insulating phases are not 
dispersed homogeneously, with the RuC>2 clusters segre- 
gated between large glassy regions of few micrometers in 
size. This segregation effect in TFRs is well known and 
it is due to the large difference in size between the fine 
conducting particles and the much coarser glassy grains 
employed in the preparation of the resistors^ Concern- 
ing the effect of the firing temperature, it is interesting 
to note that for the B series there is not much qualita- 
tive difference in the microstructure between Tf = 550° C 
[FigEfa)] and T f = 600°C [Fig0»], while for the A 
series it appears that the conducting phase is more clus- 
tered at low firing temperature [FigEf a)] than at high Tf 
[FigEfb)], where the appearance of grey regions indicate 
larger RuC>2 dispersion in the glass. 

In Figl^Jwe report the room temperature conductivity 
a measured for four different series of TFRs (see Table 
HJ as functions of the RuO-2 volume concentration x. As 
shown in Fig|5Ja), a vanishes at rather small values of x, 
as expected when the mean grain size of the conducting 
phase (40 nm and < 400 nm) is much smaller than that of 
the glass (1-5 /xm)i£2 The same data are re-plotted in the 
ln-ln plot of FigEJb) together with the corresponding fits 
to Eq.JIJ (solid lines) and the best-fit parameters <7o, x c 
and t are reported in Table [i]. As it is clearly shown, our 
conductivity data follow the power law behavior of Eq. 
with exponent t close to the universal value to ~ 2.0 for 
the Al series (t = 2.15 ± 0.06) or markedly nonuniversal 
as for the A2 case which displays t — 3.84 ± 0.06. The 
Bl and B2 series have nearly equal values of t (t ~ 3.16) 
falling in between those of the Al and A2 series. 

It is tempting to interpret the different transport be- 
haviors of the Al and A2 series by referring to the mi- 
crostructures reported in Fig[7| It appears that uni- 
versal behavior is found for the more clustered samples 
[FiglZJa)] while the nonuniversal behavior is observed 
when the conducting phase is more dispersed in the 
glass [FigUJb)]. This interpretation is coherent with the 
nonuniversality of both Bl and B2 series, which indeed 
display a large amount of RuC>2 dispersion evidenced 
by the grey regions in FigsE^a) and (b). As discussed 
in Sec[D] the microstructure has a primary role for the 
onset of nonuniversality. This is certainly true for the 
tunneling-percolation model in which the microstructure 
governs the tunneling distribution function. In this re- 
spect, the SEM images reported in Figs[7| and [H] may 
suggest that for the Al series [Fig0a)], since the R11O2 
grains are less dispersed, the tunneling distribution func- 
tion is much narrower than those of the other series. 

To measure the piezoresistive response, four resistors 




FIG. 9: (a): conductivity a as a function of RuCh volume 
concentration x for four different series of TFRs. (b) ln-ln 
plot of the same data of (a) with fits to Eq.Q shown by 
solid lines. The dashed line has slope to = 2 corresponding 
to universal behavior of transport. The prefactor po, critical 
concentration x c and transport exponent t values obtained by 
the fits are reported in Table [I] 



for each series and with equal Ru02 content were screen 
printed in a Wheatstone bridge arrangement on the top 
of alumina cantilever bars 51 mm long, 6 = 5 mm large, 
and h — 0.63 mm thick. The thermal treatment was 
the same as for the samples used to for the conductivity 
measurements. The cantilever was clamped at one end 
and different weights were applied at the opposite end. 
The resulting substrate strain e along the main cantilever 
axis can be deduced from the relation e = QMgdj (Ebti 2 ), 
where d is the distance between the resistor and the point 
of applied force, E = 332.6 GPa is the reduced AI2O3 
Young modulus, g is the gravitational acceleration and 
M is the value of the applied weight. By fixing the main 
cantilever axis parallel to the x direction, then in plain 
strain approximation the strain field transfered to the re- 
sistors is e x = e, e y = 0, and e z = — v)e, where 
v = 0.22 is the Poisson ratio of 96% A1 2 3 . Two different 
cantilevers were used for the measurements of the longi- 
tudinal and transverse piezoresistive signals obtained by 
recording the conductivity changes along the x and y di- 
rections, respectively. Then, according to Ea. (|14l) : 



Sax 
a 

5<r„ 



r„ 



1-u 



l-2v 



(32) 
(33) 



In Fig llOf a - ) and Fig llOf b') we plot the conductivity vari- 
ations along the x and y direction respectively as a func- 
tion of e for the A2 series. The RuC>2 volume fractions 
are x = 0.23, 0.154, 0.11, 0.095, and 0.085 from bottom 
to top. In the whole range of applied strains, the sig- 
nal changes linearly with e, permitting to extract from 
the slopes of the linear fits of Sai/a vs e the values 
of the longitudinal and transverse piezoresistive factors 
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TABLE I: Label legend of the various sample used in this work with fitting parameters of Eqs. 111341 



Label 


RuCh grain size 


firing temperature Tj> 


x c 


ln(<7ofim) 


t 


To 


dt/de 


Al 


400nm 


525°C 


0.0745 


11.1 ± 0.3 


2.15 ±0.06 


16.5 ±4.5 


-0.6 ± 1.2 


A2 


400nm 


600° C 


0.0670 


14.2 ± 0.2 


3.84 ±0.06 


-26.4 ±4.8 


16.2 ± 1.5 


Bl 


40nm 


550° C 


0.0626 


14.3 ±0.5 


3.17±0.16 


-45.9 ± 9 


26.1 ± 2.7 


B2 


40nm 


600° C 


0.0525 


13.7 ±0.7 


3.15 ±0.17 


-57.9 ±7.2 


33.0 ± 2.1 
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FIG. 10: (a) relative variation of conductivity along the x 
axis as a function of applied strain e in cantilever bar mea- 
surements of the A2 series for different contents x of RuC>2. 
From bottom to top: x = 0.23, 0.154, 0.11, 0.095, and 0.085. 
Solid lines are linear fits to the data, (b) the same of (a) for 
the case in which the conductivity change is measured along 
the y-axis. (c) longitudinal IV and transverse F±, piezoresis- 
tive factors obtained by applying Eqs. 1321331 to the data of 
(a) and (b), respectively. 



through Eqs . ■ The so obtained Ij and T± val- 

ues of the A2 series are plotted in Fig llOl cl as a func- 
tion of RA1O2 volume concentration x. The main feature 
displayed in Fig HTl cl is that the longitudinal and the 
transverse piezoresistive factors are not much different 
and both seem to diverge as x approaches the percolation 
threshold x c ~ 0.067. The fact that r 11 ~ T ±_ is consis- 
tent with the vanishing of the piezoresistive anisotropy 
ratio x = — F±)/Tu at the percolation threshold, as 
discussed in the previous section. 

Note instead that the transport exponent for the A2 
series is t ~ 3.84 (se Table [1J, that is much larger than 
<o = 2. Hence the divergence of T\\ and T± as x — > x c 
could well be the signature of a tunncling-percolation 
mechanism of nonuniversality. In order to investigate 
this possibility we plot in Fig lllf a') the isotropic piezore- 
sistive factor T = T|| ± 2T 1 _ for the Al, A2, Bl, and 
B2 series as a function of x. With the exception of the 
Al series which displays an almost constant piezoresis- 
tive response, the other series clearly diverge at the same 
critical concentration x c at which a goes to zero. Accord- 
ing to the tunneling-percolation theory, T should follow 
Eg. (119^1 which predicts a logarithmic divergence at the 



percolation threshold when the exponent t is nonuni- 
versal. By using the equivalence between Ea. i|17|) and 
Ea.lfTSjl our data could then be used to verify this hy- 
pothesis. This is done in Fig lllf b') where T is plotted as 
a function of ln(x — x c ) with x c values extracted from the 
conductivity data. In the entire range of concentrations, 
and for all the series Al, ...,B2, Y is rather well fitted by 
the expression: 



1 \r -f H* 



X c ) 



t = t 
t > t 



(34) 



which is Eq. I|19|) rewritten in terms of RuC>2 volume con- 
centration x. The Al series, which has DC exponent 
t = 2.15 ± 0.06 close to the universal value, has no x 
dependence of T, while the other series A2, Bl, and B2, 
characterized by nonuniversal exponents, display a loga- 
rithmic divergence of T This is in agreement 
with the expectations of the tunneling-percolaiton the- 
ory. Furthermore, as shown in Table [I] and in the inset 
of Fig lllf b). the factor To is positive for the Al series 
and negative for A2, Bl, and B2, in agreement with the 
results of the last section. For the latter series, Tq be- 
haves as r = — (1.9 ±0.5) dt/de (see inset of FiglTTTb^ in 
accord with the asymptotic relations obtained by EMA, 
Ea. (|26|) . and by our numerical calculations, Ea. (|31|l . 
The logarithmic divergence of T for the series having 
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FIG. 11: (a) Piezoresistive factor r = Tu + 2r_i_ plotted as a 
function of Ru02 volume concentration x. (b) V as a function 
of ln(a; — x c ) and fits (solid lines) to Ea. l34ll . The fit param- 
eters dt/de and To are reported in Table I and in the inset, 
together with T = — 1.9dt/de (solid line). 
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FIG. 12: (a) Piezoresistive factor V plotted as a function of 
\/(x— x c ) (symbols) with tentative fits to Eg. l35H (solid lines). 
The data related to different series have been shifted vertically 
by +30 for clarity. 



nonuniversal values of t and the corresponding negative 
values of Tq are features which can be coherently ex- 
plained by the tunncling-percolation model of nonunivcr- 
sality. However, in previous studies, this possibility was 
neglected, and the divergence of T, already reported in 
Refl37l for TFRS and in Ref [l6| for carbon-black-polymer 
composites, was attributed to a different mechanism in- 
dependent of the universality breakdown of t. This called 
into play the possibility of having nonzero derivative of 
the volume concentration x with respect to the applied 
strain when the elastic properties of the conducting and 
insulating phases are different. For the particular case 
of RuC>2-based TFRs, one finds that dx/de ~ —Ax for 
small values of the RuC>2 concentration x. It is easy to 
show that in this limit A oc 1— -B g iass/-6R U o 2 , where -B g i aS s 
and -Br,u0 2 are the bulk moduli of the glass and the con- 
ducting particles, respectively. Since -Br u o 2 — 270GPa 
and -Bgiass — 40 — 80GPa, A is expected to be different 
from zero and positive. If this reasoning held true, by 
differentiating Eq.QJ with respect to e, and by keeping t 
constant, T = —d\n(a)/de would reduce to>2£ 



r = r n 



At- 



X Xc 



X Xc 



(35) 



where we have defined K\ = Tq + At and K2 — Atx c . In 
Fig. E| we have re-plotted the T values of Fig lllf a) as a 
function of l/(x — x c ) with the same values of the criti- 
cal concentrations x c extracted from the resistivity data. 
According to Eq. I|35|) , T should follow a straight line as a 
function of l/(x — x c ) which, although being rather cor- 
rect for the A2 series, is manifestly not true for the Bl 
and B2 series. In addition, the Al series remains almost 
constant, implying that A = for this case, contrary to 
the premises of Refl37l 

Besides from the bad fit with our data, the reasoning of 
Refl3?l is based on a misunderstanding of the actual phys- 
ical meaning of x appearing in Eq.Q. In fact, x should 



be considered just as an operative estimate of the con- 
centration p of intergrain junctions with finite resistances 
present in the sample^. Current can flow from one end 
to another of the composite as long as a macroscopic 
cluster of junctions spans the entire sample. Instead of 
x, an equally valid variable describing the integrain junc- 
tion probability p would have been the concentration in 
weight of Ru02, which is manifestly independent of ap- 
plied strain (and of the e lastic p roperties of the material). 

The reasoning of RefsEH 7 is therefore non-consistent 
with the physics of percolation. However, one could ten- 
tatively argue that the applied strain actually changes 
the concentration p of junctions by breaking some of the 
bonds. This situation could be parametrized by allow- 
ing ap dependence of x so that dx/de ~ (dx/dp)(dp/de). 
Also in this case however one would end with a 1/ (x — x c ) 
divergence of T which we have seen to lead to poor fits for 
our data (Figll2[). In addition, the values of the applied 
strains in our measurements are so small (e ~ 10 ) that 
their effect is that of changing the value of the tunneling 
junction resistances without affecting their concentration 
p, so that one realistically expects that dp/de = 0. This 
is confirmed by the results of Fie llOf a) which show no 
deviation from linearity for the entire range of e values. 

In a previous publication^ we have re-analyzed the 
piezoresitive data reported in Refi3 by assuming a log- 
arithmic divergence of F rather than the l/{x — x c ) be- 
havior of Eq. H35fl ■ The agreement with Eq. i|34|) was sat- 
isfactory and, furthermore, also in this case we obtained 
negative values of Tq. However, contrary to the present 
data, all the TFRs used in Refl37l were nonuniversal, so 
it was not possible to establish the disappearance of the 
logarithmic divergence of T when t = to- 



V. CONCLUSIONS 

In this paper we have presented conductivity and 
piezoresistivity measurements in disordered Ru02-glass 
composites close to the percolation threshold. We have 
fabricated samples displaying both universal and nonuni- 
versal behavior of transport with conductivity exponents 
ranging from t ~ 2 for the universal samples up to 
t ~ 3.8 for the nonuniversal ones. The corresponding 
piezoresistive responses changed dramatically depending 
on whether the composites were universal or not. For the 
composites with t ~ 2, the piezoresistive factor T showed 
little or no dependence upon the Ru02 volume fraction 
x, whereas the nonuniversal composites displayed a log- 
arithmic divergence of T QlS X — Xc 

— > 0, where x c is the 
percolation threshold. We have interpreted the piezore- 
sistivity results as being due to a strain dependence of 
the conductivity exponent when this was nonuniversal. 
As discussed in Sec lIIII a logarithmic divergence of T 
is fully consistent with the tunncling-percolation model 
of nonuniversality proposed by Balberg a few years ago. 
According to this theory, when the tunneling distance be- 
tween adjacent conducting grains has sufficiently strong 
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fluctuations, the DC exponent acquires a dependence 
upon the mean tunneling distance a. An applied strain e 
changes a to a(l + e) which is reflected in a e-modulation 
of the exponent i, and eventually to a logarithmic diver- 
gence of r = — d\n{a)/de at x c . By studying an effec- 
tive medium approximation of the tunncling-percolation 
model, and by extensive Monte Carlo calculations, we 
have shown that when T diverges, the x-independent con- 
tribution Tq of r becomes negative, in agreement with 
what we observed in the experiments. 

In view of such agreement between theory and exper- 
iments, and given the fact that an alternative explana- 
tion based on a strain dependence of the integrain junc- 
tion concentration is non-physical and leads to poor fits 
to our data, we conclude that the origin of nonuniver- 
sality in Ru02-glass composites is most probably due 
to a tunneling-percolation mechanism of nonuniversality. 
This conclusion is also coherent with the observed corre- 
lation between the onset of nonuniversality and the mi- 
crostructure of our samples, which showed an highly clus- 
tered arrangement of the conducting phase when t = t a 
or a more dispersed configuration when t > to- Although 
being only qualitative, this picture suggests a possible 
route for more quantitative analysis on the interplay be- 
tween criticality and microstructure. 

The tunncling-percolation mechanism of universality 
breakdown could apply also to other materials for which 
transport is governed by tunneling such as carbon-black- 
polymer composites, and experiments on their piezoresis- 
tive response could confirm such conjecture. Some ear- 



lier data showing diverging piezoresistivity response at 
the conductor-insulator transition already exist jiLZi but 
their interpretation is not straightforward due to non- 
linear conductivity variations as a function of strain or 
pressure, or even to hysteresys effects. 

An interesting issue we have not addressed in the 
present work is the possible effect of temperature T 
on the piezoresistivity of disordered composites. As 
shown in Refl&iL at fixed concentration of the conduct- 
ing phase, carbon polyvinylchloride composites display 
a rather strong enhancement of T at low T. This was 
interpreted in terms of thermally activated voltage fluc- 
tuations across the tunneling barriers^ A qualitatively 
similar enhancement of T as the temperature drops is 
expected also in models based on variable-range hopping 
mechanism of transport »S It should be however pointed 
out that in these works the problem of connectivity is not 
included so that the tunneling current flows on a net- 
work without percolation characteristics. Nevertheless, 
the percolation picture (with its corresponding critical 
exponents) seems to survive well also in the low tempera- 
ture tunneling regime^. Hence, a tunneling-percolation 
theory of piezoresistivity at low temperatures should be 
formulated by considering percolation networks where 
the simple tunneling process, Eq.@, are combined with 
additional terms describing, e. g., grain charging effects 
and/or Coulomb interactions. 

This work was partially supported by TOPNANO 21 
(project NAMESA, No. 5557.2). 
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